Subharmonic Shapiro steps and assisted tunneling in superconducting point contacts 
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We analyze the current in a superconducting point contact of arbitrary transmission in the presence 
of a microwave radiation. The interplay between the ac Josephson current and the microwave signal 
gives rise to Shapiro steps at voltages V = {m/n)hujr/2e, where n,m are integer numbers and Ur 
is the radiation frequency. The subharmonic steps (n 7^ 1) are a consequence of multiple Andreev 
reflections (MAR) and provide a signature of the peculiar ac Josephson effect at high transmission. 
Moreover, the dc current exhibits a rich subgap structure due to phot on- assisted MARs. 



Introduction. — Our understanding of the electronic 
transport through superconducting nanostructures has 
experienced a notable development in last few years 0. 
Partly, this has been due to the appearance on scene of 
the metallic atomic-size contacts, which can be produced 
by means of scanning tunneling microscope and break- 
junction techniques ||-|^. These nanowires have turned 
out to be ideal systems to test the modern transport the- 
ories in mesoscopic superconductors. Thus, for instance 
Scheer and coworkers ||] found a quantitative agreement 
between the measurements of the current- voltage charac- 
teristics of different atomic contacts and the predictions 
of the theory for a single-channel superconducting con- 
tact These experiments not only helped to clarify 
the structure of the subgap current in superconducting 
contacts, but also showed that the set of the transmis- 
sion coefficients in an atomic-size contact is amenable to 
measurement. This possibility has recently allowed a set 
of experiments that confirm the theoretical predictions 
for transport properties like supercurrent [Q and noise 
From these combined theoretical and experimental 
efforts a coherent picture of transport in superconducting 
point contacts has emerged with multiple Andreev reflec- 
tions (MAR) as a central concept. However, in spite 
of these recent successes, one of the most remarkable pre- 
dictions of MAR theory remains to be confirmed, namely 
the ac Josephson effect. The theory says that in a con- 
stant voltage biased superconducting point contact, the 
time-dependent current is given by I{t) = -^"S™'^"*- 
This means that the occurrence of MARs gives rise to 
the appearance of alternating currents that oscillate not 
only with the Josephson frequency luq = 'ieV/Ti, V being 
the voltage, as in the case of tunnel junctions, but also 
with all its harmonics. So far there is no experimental 
evidence of the existence of such components. 

In this paper, we present a theoretical analysis of the 
current in a superconducting point contact under a mi- 
crowave radiation. We show that the interplay between 
the ac Josephson current components and a microwave 
signal leads to the appearance of Shapiro steps at voltages 
V = (rn/n)hLiJr/'2.e, where n, m are integer numbers and 
ijjr is the frequency of the radiation. This means that in 



addition to the usual steps {n = 1) found in tunnel junc- 
tions there also appear subharmonic Shapiro steps 
{n ^ I), which constitute an unambiguous signature of 
the ac Josephson effect in these contacts. Moreover, we 
also find that the dc background current, in which the 
Shapiro steps are superimposed, exhibits a rich subgap 
structure, which can be understood in terms of photon- 
assisted MARs and provides a natural explanation of ex- 
perimental findings in the early seventies |l0[ . 

Theoretical model. — Our goal is to calculate the cur- 
rent in a voltage biased superconducting quantum point 
contact (SQPC) [0 in the presence of a monochromatic 
radiation of frequency ujr. We assume that the external 
radiation produces an effective time-dependent voltage 
V{t) = V + Vac sin ujrt. Our task is to extend the MAR 
theory to the case of such a time-dependent voltage, for 
which the so-called Hamiltonian approach is a conve- 
nient starting point. For the voltage range eV A one 
can neglect the energy dependence of the transmission 
coefficients and all transport properties can be expressed 
as a superposition of independent channel contributions. 
Thus, the problem reduces to the analysis of a single 
channel contact, which can be described by means of the 
following tight-binding-like Hamiltonian 



H = Hl + Hr 



E{-1 



}, (1) 



where Hl,b. are the BCS Hamiltonians for the isolated 
electrodes. In the coupling term L and R stand for the 
outermost sites of each electrode, and w is a hopping 
parameter coupling these sites. This parameter deter- 
mines the normal transmission coefficient of this model 
T, which adopts the form T = A{v/W fl [l + [v/Wf] ^ 
where W — Xji^pp, with pp being the electrodes density 
of states at the Fermi energy 

In this model the current evaluated at the interface 
between the two electrodes adopts the form 



= i E {^(^Lwci^^w) - «*(cLwcL.(i))} . (2) 
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The non-equilibrium expectation values in Eq. (g) can be 
expressed in terms of the Keldysh Green functions Gfj~ , 
which in the 2x2 Nambu representation read 



Thus, the current can be now written as 



m = Tr 



T3 



{v{t)G^Z{t,t)-v\t)Gt],it,t))] , (4) 



where is the corresponding Pauli matrix, Tr denotes 
the trace in Nambu space and v is the hopping that in 
the Nambu matrix representation is written as 
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(5) 



Fourier components, T^^^ = T{e + neV + kfiujr, e + meV- 
IfuVr), in the following way 
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(7) 



At this point, the calculation of the current has been 
reduced to determination of the Fourier components of 
the T-matrix. In the case of a symmetric contact con- 
sidered here, one can show that the dc current can be 

expressed only in terms of 
following set of linear algebraic equations 



^Lfi.iO' which fulfill the 



Here, (/)(t) — (pQ + uj^t + 2acos Wr-t is the time-dependent 
superconducting phase difference. The constant a = 
eVacI {hiOr) measures the strength of the coupling to the 
electromagnetic field, and is proportional to the square 
root of the radiation power. 

In order to determine the Green functions we follow 
a perturbative scheme and treat the coupling term in 
Hamiltonian (1) as a perturbation. The unperturbed 
Green functions, g, correspond to the uncoupled elec- 
trodes in equilibrium. Thus, the retarded and ad- 
vanced components adopt the BCS form: ^'"'''(e) = 
g'^'^ie)! + /''■°(e)-ri, w here g'~''^(e) = -(e'''''/A)/(e) = 
-e'^'^/W^J/S? - (e''''^)2, where e"'" = e±iri, with 77 = 0+. 
Following Ref. § we express the current in terms of a T- 
matrix, rather than in terms of the Green functions. The 
T-matrix associated to the time-dependent perturbation 
of Eq. (5) is defined as T'"''^ = v + vo'g^'°-oT^'^, where the 
o product is a shorthand for integration over intermedi- 
ate time arguments. As shown in Ref. |^], the current in 
terms of the T-matrix components reads 

m = I Tr [f3 (tIj, o g+R- o nL ogl-glo ° 9^ o f^^ 
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where the different matrix coefficients adopt the following 
form in terms of the unperturbed Green functions 

\ 9i-iJi 9i-i9i . 



^i,i+2 — ^ ' 







V 9i-2 Ji-2 

J 



+9r o t^bl ° ° tIr -T^RLogt-o o g%) \ , (6) 

In order to solve the T-matrix integral equa- 
tion it is convenient to Fourier transform with 
respect to the temporal arguments, T{t,t') = 



Iq is a back- 



(l/27r)/de Jde' e' 



r(e, e'). Due to time depen- 



dence of the coupling element (see Eq. (5)), one can show 
that ^(e, e') admits the following solution: T{e, e') = 
J2n m'^i^'^ + "-^^ + m^^r)5{e — e' + neV + mfiLUr)- 
Thus, one can finally write down the current as I{t) — 
m -^rT ^^P I* (^"/"o + nujot + mujrt)], where the current 
amplitudes /™ can be expressed in terms of the T-matrix 



where we have used the shorthand notation = g"'{^ + 
ieV + khojr) and Jnio:) is the Bessel function of order 
n. In some limits one can find an analytical solution of 
these systems, but in general a numerical calculation is 
needed. 

Results and discussions. — Let us concentrate in the 
dc current, Idc- This current is the sum of two contri- 
butions: Idc = Ib + I Shapiro, wherc Ib 
ground current and Ishapiro = 
is the Shapiro steps contribution at discrete voltages 
= {m/n)hujr/2e. Notice that several ac current am- 
plitudes can give a dc contribution at the same voltage. 
Notice also that the Shapiro step contribution depends on 
the average value of the phase, (pQ. We shall concentrate 
in the height of the Shapiro steps, which will be denoted 
as 5,™. Let us remark that in the tunneling regime we 
recover the well-known results for both the background 
current and Shapiro step heights |T2[ . 

In order to illustrate the general results, in Fig. 1 we 
show the dc current, background current plus Shapiro 
steps, for different values of a and a frequency ujr = 0.5A. 
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We can see the two main features that wiU be the sub- 
ject of the rest of the paper: (i) the subharmonic Shapiro 
steps 5™, with n ^ 1, are clearly visible at high trans- 
missions, and (ii) the background current exhibits a sub- 
harmonic gap structure at voltages eV = {2A + kh(jJr)/n, 
with n, k integers, which is specially pronounced at low 
transmissions llSl. 




0.5 1 1.5 2 2.5 0.125 0.25 0.375 0.5 



eV/A eV/A 

FIG. 1. Zero temperature dc current, Idc, as a function 
of voltage for a frequency Ur = 0.5A and several values of 
a. The different curves in each panel correspond to different 
transmissions, as indicated in panel (a). In panels (b) and (c) 
the curves have been vertically displaced. Panel (d) shows 
in detail the curve T = 0.95 of panel (c). The current is 
normalized by the normal conductance Gn = (2e^//i)T. 

Let us start by analyzing the background current. In 
Fig. 2(c-d) we show the background current for two dif- 
ferent frequencies at a moderate power, a — 1.0. The cur- 
rent in absence of radiation is also shown for comparison. 
As mentioned above, the most prominent feature in the 
background current is the appearance of a pronounced 
subgap structure at voltages eV = (2A + khujr)/n. This 
structure is specially clear at low transmissions (see Fig. 
2d) and progressively disappears as the transparency is 
increased. Indeed, this peculiar subharmonic gap struc- 
ture was already observed in several experiments in the 
early seventies in point contacts and thin-film micro- 
bridges |l^. At that time no consistent explanation 
was given, but it is clear that this structure can be ex- 
plained in terms of photon-assisted MARs. A step at 
eV — (2A -|- kujr)/n is simply due to the opening of a 
MAR of order n in which k photons in total are ab- 
sorbed (fc negative) or emitted (fc positive). This is il- 
lustrated in the upper panels of Fig. 2. In order to 
understand how this subharmonic structure evolves with 
the rf power, one can do a systematic perturbative ex- 
pansion in the transmission. This analysis tells us that 
at low transparency the height of a current jump at 



eV = (2 A -|- ku!r)/n is proportional to J'^{na), which 
is valid as long as huJr <?C 2A/n. This results coincides 
with the phenomenological functional form that was used 
to fit the experiments by Soerensen et al. 
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FIG. 2. (a) Representation of a MAR of order 3 in which 
k photons are absorbed. This process has a threshold volt- 
age eVth ~ (2A — h\k\u}r)/-i, and its probability amplitude 
is proportional to Jk[ot)- (b) A 3-order MAR mediated by 
the emission of k photons, which contributes to the subgap 
structure at eVth ~ {2A + h\k\LJr) /3. (c) Background current 
as a function of voltage for cur = 0.5A and different transmis- 
sions, (d) The same as in (c) for ujr = O.IA. The inset shows 
a blow up around eV — A. The dotted lines in (c) and (d) 
correspond to the current in absence of radiation. 

Let us now discuss the Shapiro steps. In this case 
the most important aspect is the existence of subhar- 
monic steps absent in tunnel junctions. These steps arise 
from the phase locking between the harmonics of the 
Josephson frequency and the harmonics of the ac radi- 
ation. Early experiments on the ac Josephson effect in 
weak links observed subharmonic steps in the I-V curves 
[Q. More recently, there have been reported observa- 
tions of non-integer Shapiro steps in high-Tc contacts 
|]l5| , S-semiconductor-S junctions [|6| and diffusive S- 
N-S systems 0. Although the Shapiro steps can be 
understood as a simple consequence of a non-sinusoidal 
current-phase relation, the present approach goes beyond 
a simple "adiabatic" approximation and provides the first 
microscopic theory of Shapiro steps in contacts of arbi- 
trary transmission. The adiabatic approximation, which 
introduces the time-dependence into the zero bias super- 
current through the Josephson relation, gives rise to the 
well known Bessel-function-like behavior of the steps and 
gives a good descprition of the tunnel regime |l2j . How- 
ever, as we show below, such a simple approach fails in 
the description of a highly transmissive contact. 

As a rule of thumb, a Shapiro step is visible when 
the corresponding ac Josephson component, /„, in ab- 
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sence of radiation gives a significant contribution. In 
particular, this means high transmissions (see Figs. 3- 
4 in Ref. 0). One can show that the leading order in 
transmission of a Shapiro step goes like ~ T", which 
is a consequence of the fact that /„ ^ T" , and the reason 
for the absence of the n ^ 1 steps in low transmissive 
contacts. However, near perfect transmission the sub- 
harmonic steps can be even higher than the integer ones. 
This behavior is illustrated in Fig. ^, where we show the 
Shapiro steps S'^ as a function of the transmission for two 
different frequencies. 
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FIG. 3. Shapiro steps S\ versus transmission for a — 0.25. 



Fig. ^ shows the power dependence of the Shapiro 
steps for a frequency uj^ — 0.5A. Notice that this 
dependence is rather complicated for both integer and 
subharmonic steps, and clearly deviates from the usual 
Bessel function behavior. This is due to the frequency- 
dependence of the Josephson components, which is spe- 
cially pronounced at high transmissions. Neglecting this 
dependence, i.e. within an adiabatic approximation, one 
would get that S"™ evolves as |Jm(2na)|. However, as 
shown in Fig. as the transmission increases the va- 
lidity of this approximation is restricted to a ^ 1. Notice 
also the complex oscillation pattern at high transmissions 
(see T = 0.8 curves in Fig which is due to the fact 
that several ac components give a significant contribution 
to the same Shapiro step. 

In summary, we have presented a theoretical analysis of 
the dc current in a superconducting point contact in the 
presence of a microwave radiation. We have shown that 
the microscopic theory of coherent multiple Andreev re- 
flections provides an unified description of Shapiro steps 
and assisted tunneling, explaining in a natural way the 
observations of subharmonic steps fl^-lr^ and the pecu- 
liar subharmonic gap structure under a microwave radia- 
tion [0 . Let us finally remark that the results presented 
in this work are amenable to a quantitative experimental 
test using atomic-size contacts 
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FIG. 4. (a-c) Shapiro steps S\ (n — 1,2,3) as a function 
of a for u)r = 0.5A. The dotted lines in panel (b) correspond 
to the adiabatic approximation: ~ |Ji(4a)|. 
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